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The recently discovered 150-cavity (formed by loop residues 147-152, N2 numbering) adjacent to the 
enzymatic active site of group 1 influenza A neuraminidase (NA) has introduced a novel target for the design 
of next-generation NA inhibitors. However, only group 1 NAs, with the exception of the 2009 pandemic 
H1N1 NA, possess a 150-cavity, and no 150-cavity has been observed in group 2 NAs. The role of the 
150-cavity played in enzymatic activity and inhibitor binding is not well understood. Here, we demonstrate 
for the first time that oseltamivir carboxylate can induce opening of the rigid closed N2 1 50-loop and 
provide a novel mechanism for 1 50-loop movement using molecular dynamics simulations. Our results 
provide the structural and biophysical basis of the open form of 1 50-loop and illustrates that the inherent 
flexibility and the ligand induced flexibility of the 1 50-loop should be taken into consideration for future 
drug design. 

I nfluenza virus causes epidemics and pandemics, which severely impair public health 15 . Two of the major 
I pandemics of the last century were caused by N2 containing influenza viruses: H2N2 (Asian flu) and H3N2 
I (Hong Kong flu) 1 . Recent outbreaks of H3N2 and H1N2 swine influenza viruses-caused human febrile 
respiratory illness in the United States highlight the importance of this subtype in influenza epidemics 6 8 . 
Besides H1N1, the H3N2 subtype is the major causative agent of severe epidemics and is critical for vaccine 
development 9 . Therefore a comprehensive understanding of N2 containing viruses is important for preparedness 
against highly transmissible influenza viruses. 

Hemagglutinin (HA) and neuraminidase (NA) are the two major surface glycoproteins responsible for ini- 
tiating influenza virus infection 1013 and virus release 1416 , respectively. HA and NA of influenza A viruses are 
divided into subtypes based upon their distinct antigenic properties: seventeen for HA (H1-H17) and ten for NA 
(N1-N10) 17 " 19 . Among the influenza A viruses, only Nl and N2 have been found in human isolates responsible for 
pandemics and recurrent annual epidemics. With the exception of N10, recently identified in a bat influenza A 
virus genome, the nine NA alleles are classified into two groups according to phylogenetic analysis and structure. 
Group 1 NA comprises Nl, N4, N5 and N8, whereas group 2 comprises N2, N3, N6, N7 and N9 20 . 3-D structures 
reveal the distinct conformations of the areas adjacent to the enzymatic active site between group 1 and group 2 
members though the active site structures are virtually identical among all the NAs 20 . For typical Nl subtypes (but 
not for the 2009 H1N1 pandemic Nl), the crystal structures reveal a 150-loop (formed by amino acids 147-152, 
N2 numbering) that adopts an open conformation forming an additional 150-cavity adjacent to the active site. 
Previously, no group 2 NAs have been crystallographically shown to have a 150-loop in an open conformation, 
although all-atom molecular dynamics simulations indicate that N2 may be able to adopt such a configuration in 
solution 21 . The recently discovered 150-cavity is currently being explored as a novel target for group 1 specific 
influenza NA inhibitors 20,22 . 

Detailed structural analyses demonstrate that residue 147 plays an essential role in the conformation of the 
150-loop. Recently, our group has shown that the N5 structure contains an extended 150-cavity resulting from the 
unique residue N147 23 . In the solved N2 structures, a salt bridge between D147 and H150 contributes to a rigid 
closed 150-loop 24 . Previous molecular dynamics simulations show that the D147-H150 salt bridge greatly 
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stabilizes the closed 150-loop conformation and that loss of this salt 
bridge reduces the rigidity of the 150-loop 21 . D147 is prevalent in N2, 
but rarely found in other types of influenza NA 23 . Compared with N2 
containing D147, the flexibility of the 150-loop of 2009 H1N1 NA 
(G147) is much higher 21 , although the structure of 2009 H1N1 NA 
also presents a deficient 150-cavity in its crystal structure 25 . 

Since NA plays an essential role in the release of new influenza 
virions from host cells, inhibition of NA compromises the ability of 
progeny virions to spread to uninfected cells. Oseltamivir (Tamiflu) 
and zanamivir (Relenza) are two commercially available NA-tar- 
geted competitive inhibitors, which act against both group 1 and 
group 2 enzymes as well as influenza B NA 16 . The open 150-loop 
of group 1 NAs has been found to adopt a closed conformation upon 
binding of zanamivir, however for oseltamivir carboxylate this effect 
depends on both inhibitor concentration and soaking time depend- 
ent 20 . The crystal structures of typical group 1 NAs in complex with 
oseltamivir carboxylate display two 150-loop conformations, 
indicating a two-step process of oseltamivir carboxylate binding. 
Molecular dynamics simulations of the free and oseltamivir carbox- 
ylate-bound forms of tetrameric Nl suggest a rapid loop switching 
motion, which demonstrates the flexibility of the 150-loop 26 . In 
contrast, in the crystal structures of typical group 2 NA-inhibitor 
bound complexes, the 150-loops always adopt closed conforma- 
y on 2o,27-2o Also, it is noteworthy that zanamivir always induces the 
closed conformation 20,27 , whether the target is a group- 1 or group-2 
NA. Moreover, previous reports showed that R152K mutant in flu B 
lead to zanamivir and oseltamivir resistance 30,31 , which suggests that 
150-loop may also play a role in inhibitor binding. 

Here we report a novel half open 150-loop in the crystal structure 
of a typical group 2 NA for the first time. This suggests that inhibitors 



targeting the 150-cavity may also effective target group 2 influenza 
NAs. Furthermore, we demonstrate with molecular dynamics simu- 
lations that the protonation state of inhibitors and N2 HI 50 play an 
important role on the movement of the 150-loop. Comparison of the 
differences between oseltamivir carboxylate and zanamivir binding 
to N2 illustrates the importance of the inhibitor C4 group (zanamivir 
numbering), providing considerable value for influenza inhibitor 
design. 

Results 

Oseltamivir carboxylate induced opening of the rigid N2 150- 
loop: the first open 150-cavity structure of group 2 NA. Se- 
quence alignment shows that D147 is commonly found in N2 21,23 . 
D147 can form a rigid salt bridge with H150, which plays an impor- 
tant role in holding the N2 150-loop in a closed conformation (PDB 
1NN2) 24 . Thus, all the N2 crystal structures in the PDB database 
contain a rigid closed 150-loop. The free N2 observed in our 
experiment also presents a closed 150-loop (hereafter referred to as 
"Free N2") (Fig. la). 

However, upon soaking with oseltamivir carboxylate, we found 
that the rigid N2 150-loop can actually be opened (Fig. la). Our 
crystal structure data indicates that the conformation of the N2 
150-loop is dependent on the concentration of oseltamivir carbox- 
ylate (Fig. la). Specifically, N2 presents an open 150-cavity when 
soaking crystals in 20 mM oseltamivir carboxylate for 1 hour (here- 
after referred to as "N2-oseltamivir open form"), while it shows the 
same conformation as in the free N2 when soaking into 40 uM 
oseltamivir carboxylate for 30 min (hereafter referred to as "N2- 
oseltamivir closed form") (Fig. la). 
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Figure 1 | Comparison of oseltamivir carboxylate induced 150-loop conformational change in N2, Nl and N8. (a) The active sites of free N2 
(yellow), Nl (cyan), N8 (lightblue) and oseltamivir carboxylate-bound complex in different soaking conditions are displayed in surface representation. 
Free N2 (yellow), N2-40 uM-30 min (orange), Nl-20 mM-3 days (green) and N8-500 uM-30 min (pink) have no 150-cavity, however free Nl, N8, 
N2-20 mM-1 hour (limon), Nl-20 mM-30 min (marine) and N8-20 |rM-150 min (magenta) contain a 150-cavity. Oseltamivir carboxylate is shown in 
wheat stick, (b) Superposition of the 150-loop of free Nl, N2, N8 and Nl, N2, N8 complexed with oseltamivir carboxylate under different soaking 
conditions. The 150-loops are shown in cartoon representation, and colors are corresponding to Fig. la. 
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An influence of oseltamivir carboxylate concentration and soak- 
ing- time on the 150-loop conformation has also been observed in 
both Nl and N8 (Fig. la), indicating that oseltamivir carboxylate has 
great influence on the conformation of the 150-loop. When using 
typical group 1 NA Nl (PDB 2HTY) and N8 (PDB 2HT5) 20 and 
typical group 2 N2 as open and closed references, it is clear that 
the 150-loop of N2 -oseltamivir open form adopts a half-open con- 
formation (Fig. lb), which is distinct from other types of NA. 
Superposition of the free N2 Cot atoms with those of N2- oseltamivir 
open form and N2- oseltamivir closed form reveals root-mean- 
square deviations (RMSD) of 3.0 A and 0.1 A in the 150-loop, 
respectively (only Cot atoms were used to calculate RMSD). The main 
chain of the 150-loop in N2-oseltamivir closed form is the same as 
free N2. The 150-loop conformation in the oseltamivir carboxylate 
complexes with Nl and N8 is similar to free N2 with a RMSD of 0.5 A 
and 0.6 A in the 150-loop, respectively. 

Intermolecular interaction of the tetrameric NA plays an essential 
role in stabilizing the open N2 150-cavity. In the N2-oseltamivir 
open form crystal structure the novel half-open form of the 150- 
cavity appears to be quite stable. The side chain of D 147 is fixed by 
the side chain of R107 in the neighboring NA molecule and a water 
molecule forms hydrogen bonds with D147 and HI 50 with distances 
of 2.75 A and 2.79 A, respectively (Fig. 2a). Additionally, the 150- 



loop turns into a 3 10 helix, further stabilizing an open 150-loop 
conformation (Fig. 2b). To investigate the importance of the D147- 
R107 intermolecular interaction in maintaining the open conforma- 
tion of the N2 150-loop, D147 was replaced by glycine. Although the 
crystal of N2-G147 was soaked in the same condition as N2- 
oseltamivir open form, the 150-loop of the N2-G147-oseltamivir 
carboxylate complex still adopts a closed conformation identical to 
that of free N2 with D147 (Fig. 2c). Specifically, the glycine does not 
have a side chain, so it cannot form hydrogen bond with R107 in 
neighboring molecule to further stabilize the 150-loop in an open 
conformation. Still, the side chain of N2-G147 H150 is flexible and 
adopts two conformations, in one of these conformations the H150 
ND1 forms hydrogen bond with the oxygen atom of the G147 main 
chain, stabilizing a closed 150-loop conformation (Fig. 2c). Addi- 
tionally, the interaction between the side chain of D151 and the 
amino group of oseltamivir carboxylate is another factor that stabi- 
lizes the 150-loop in a closed conformation. In the N2-oseltamivir 
closed form structure, the amino group of oseltamivir carboxylate 
forms a hydrogen bond with the side chain of D151, at a distance of 
3.26 A (Fig. 2d). While in the N2-oseltamivir-open structure, there is 
no interaction between the OD1 on the side chain of D151 and the 
amino group of oseltamivir carboxylate, with a distance of 8.23 A 
(Fig. 2d). 
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Figure 2 | Structural basis of an open 150-loop in N2 neuraminidase, (a) Overview of the open form of 150-loop in the N2-oseltamivr complex structure 
with an emphasis on the intermolecular interaction indicated with a violetpurple box. The residues of N2 -wild type are shown in limon, while the residues 
of N2-G147 are shown in wheat. The water molecule is shown in red sphere, (b) Superposition of the 150-loop conformation between free N2 (yellow), 
N2-oseltamivir closed form (orange) and N2-oseltamivir open form (limon). Compare the interaction between D147 and H150. 
(c) Superposition of 150-loop conformation between N2-oseltamivir closed form (orange) and N2-G147 after a 1 hour soak with 20 mM oseltamivir 
carboxylate (wheat). Compare the interaction between Residue 147 and 150. (d) Superposition of 150-loop conformation between N2-oseltamivir closed 
form (orange) and N2-oseltamivir open form (limon). Compare the interaction between D151 and amino group of oseltamivir carboxylate. The side 
chain of the residues is displayed in stick representation with hydrogen bonds and salt bridges indicated by dotted lines. The interaction with a distance 
greater than 3.5 A is shown as a broken line. The shared residues are highlighted in black and labeled. 
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Table 1 Information of five 
Systems 


simulation systems 
Ligand 


Charge in G4 group (zanamivir numbering) of inhibitor 


Protonation state of H 1 50 


OST-NH3-H150-Pos* 


Oseltamivir carboxylate 


Positive 


Positive 


OST-NH2-H150-Pos 


Oseltamivir carboxylate 


Neutral 


Positive 


OST-NH2-H150-Neu 


Oseltamivir carboxylate 


Neutral 


Neutral 


FREE-H 1 50-Neu 




Neutral 


ZMR-H1 50-Pos** 


Zanamivir 


Positive 


Positive 


*OST refers to oseltamivir carboxylate. 
**ZMR refers to zanamivir. 



Mechanism for oseltamivir carboxylate induced opening of the N2 
150-loop. The crystal structures of N2 in complex with oseltamivir 
carboxylate and zanamivir were determined at pH 9.0, while the pKa 
values of the oseltamivr amino group and zanamivir guanidino 
group are 7.75 and 13 32 , respectively. Thus, the amino group of osel- 
tamivir carboxylate should be neutral approximately 94.7% of the 
time while the guanidino moiety in zanamivir should be positively 
charged almost 99.9% of the time according to the Henderson- 
Hasselbach equation. The pKa of HI 50 in the closed 150-loop con- 
formation is 1 pK unit higher than that in the open conformation as 
predicted by H++ 33 . This indicates that HI 50 is more likely to be 
protonated in the closed conformation compared to the open confor- 
mation. To explore the mechanism of how inhibitors influence con- 
formational changes of the 150-loop we performed three separated 
200 ns molecular dynamics simulations (MDs) for tetrameric N2 
with different protonation states of HI 50 and the amino group of 
oseltamivir carboxylate. One 200 ns MDs with zanamivir and one 
200 ns MDs without any ligand were also performed (Table 1). 

In the OST-NH3-H150-Pos system, the C4 amino group of osel- 
tamivir carboxylate and HI 50 are all positively charged. The 150- 
loop remains closed for 98.1% of the conformations sampled in the 
simulation (Table 2, Supplementary Fig. SI). The positively charged 
oseltamivir carboxylate amino group forms a hydrogen bond with 
the D151 side chain, which stabilizes the closed 150-loop conforma- 
tion. The salt bridge between D147 and HI 50 contributes greatly to a 
closed 150-loop (Supplementary Fig. SI). 

In the OST-NH2-H1 50-Pos system, the oseltamivir carboxylate 
C4 amino group is neutral and HI 50 is positively charged. During the 
simulation, an open 150-cavity can be observed in all four chains and 
40.8% of the configurations present a 150-cavity volume larger than 
40 A 3 (Table 2). Specifically, in chain A and chain D, the 150-loop 
started from a closed form, opened for a while, resulting in the 
formation of a 150-cavity, and then closed again, leading to absence 
of the 150-cavity. In the other two chains the 150-loop also opened 
(Supplementary Fig. S2). The open and closed conformation of the 
150-loop are directly related to the interaction between D147 and 
H150 (Supplementary Fig. S2). Compared with the OST-NH3- 
H150-Pos system, the deprotonated amino group of oseltamivir car- 
boxylate exerts a significant effect on 150-loop opening. The inter- 
action between the amino group of oseltamivir carboxylate and 
the side chain of D151 is almost completely lost (Fig. 3a) which 
results in increasing flexibility of the 150-loop. The open and closed 



Table 2 | Population 


analysis of 1 50-cavity 




System 


Volume 


rate (%) 


0-40 


>40 


OST-NH3-H1 50-Pos 


98.1 


1.9 


OST-NH2-H1 50-Pos 


59.2 


40.8 


OST-NH2-H1 50-Neu 


39.9 


60.1 


FREE-H 1 50-Neu 


72 


28 


ZMR-H1 50-Pos 


97.8 


2.2 



conformation of the 150-loop and environmental change of H150. 
Histidine is sensitive to environmental changes, so when oseltamivir 
carboxylate influences the rigid closed 150-loop structure, H150 may 
lose a proton which will reduce the interaction between D147 and 
H150, and further lead to the 150-loop open. 

To validate the above hypothesis, we performed a MD simulation 
for the OST-NH2-H1 50-Neu system. In this system, both the amino 
group of oseltamivir carboxylate and the HI 50 imidazole side chain 
are in neutral states. A statistical analysis shows that the 150-cavity 
with a volume larger than 40 A 3 occurs during 60.1% of the simu- 
lation (Table 2), which is 19.3% and 58.2% larger than the OST-NH2- 
H150-Pos system and OST-NH3-H150-Pos system, respectively. 
The 150-loop opening is accompanied by the loss of the D147- 
H150 interaction in all four chains, although different behaviors 
can be observed in different chains and at specific times. In chain 

A, the 1 50-cavity volume decreased due to the orientation of R430. In 
chain B, the side chain of D147 failed to interact with R107 in the 
neighboring molecule, and then D147 and H150 formed weak 
hydrogen bonds again leading to the reformation of a closed 150- 
loop conformation (Supplementary Fig. S3). In chain D, the 150-loop 
opening process was observed in a similar manner observed in crystal 
soaking experiments (Fig. 3b). After 40 ns, the interaction between 
the oseltamivir carboxylate amino group and D151 was lost (Fig. 3b) 
enabling the 150-loop to be more flexible. The interaction between 
D147 and H150 is reduced as a result of deprotonation of H150, 
thereby increasing the flexibility of D147. Then the side chain of 
R107 in the neighboring molecule has the opportunity to form salt 
bridge with D147 enabling the 150-loop to form a stable open con- 
formation. 

To validate the effect of the protonation state of HI 50 on the 
flexibility of the 150-loop, the FREE-H150-Neu system with H150 
in the neutral state was set up. Instantaneous opening of the 150- 
cavity can be observed in all four of the chains. Specifically, in chain 

B, C and D, the 150-cavity opened for a short time and then closed 
again. For chain A, the cavity remains open for a longer time 
(Supplementary Fig. S4). 72% of the conformations have a 150-cavity 
smaller than 40 A 3 , which is much larger than the closed conforma- 
tion population in OST-NH2-H150-Neu (39.9%) (Table 2). This 
indicates that without oseltamivir carboxylate, N2 prefers a closed 
150-loop conformation. Although the interaction between D147 and 
H150 becomes weaker as the deprotonation of H150, hydrogen 
bonds such as R156-D151, Q136-H150 and T148-T439 also contrib- 
ute to a closed 150-loop conformation (Fig. 3a). 

The guanidino group of zanamivir is virtually always positively 
charged due to its high pKa. Thus, in the ZMR-H1 50-POS system, 
the guanidino group and HI 50 are both set as positively charged. 
During the simulation, 97.8% conformations sampled in this traject- 
ory has a 1 50-cavity smaller than 40 A 3 . The 1 50-loop remains closed 
almost throughout the entire simulation (Supplementary Fig. S5), 
and has strong interactions with zanamivir. Compared with oselta- 
mivir carboxylate, the hydrogen bond occupancy between zanamivir 
and D151 is stronger (Fig. 3a), which indicates that zanamivir has a 
higher avidity with the 150-loop. Moreover, the N2-inhibitor com- 
plex structure clearly shows that the interaction between the 
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Figure 3 | Driving force for the opening of the N2 150-loop. (a) Hydrogen bond occupancy related to the 150-loop in the five simulation systems. Grey: 
FREE-H150-Neu; green: ZMR-H150-Pos; yellow: OST-NH3-H150-Pos; pink: OST-NH2-H150-Pos; red: OST-NH2-H150-Neu. (b) The opening 
process of the 150-loop in OST-NH2-H150-Neu trajectory (Chain D) represents the similar conformational change to the crystal structure. Upper panel: 
150-cavity volume; down: the distance between alpha carbon of D151 and nitrogen atom of Oseltamivir carboxylate amino group (black), the nearest 
distance between OD1 or OD2 of D147 and ND1 of H150 (green), and the nearest distance between OD1 or OD2 of D147 and NE atom of R107 in the 
neighboring molecule (red), (c) Comparison of the interactions between N2 residues and amino group of oseltamivir carboxylate and guanidino group of 
zanamivir (within 3 A). 



zanamivir guanidino group and neighboring residues (less than 
3. 5 A) is stronger than that of the oseltamivir carboxylate amino 
group (Fig. 3c). Specifically, the zanamivir guanidino can bind to 
four residues (El 19, D151, W178 and E227), while the amino group 
of oseltamivir carboxylate can only interact with two residues (El 19 
and D151). The interactions between zanamivir and the main chain 
of the 150-loop also help to maintain the closed 150-loop conforma- 
tion. The structure and simulation data suggest that the proper size 
and charge of the C4 moiety plays a critical role in determining 150- 
loop dynamics upon ligand binding. 

Discussion 

The influenza NA 150-cavity was first observed in group 1 NAs, and 
several group 1 specific inhibitors were subsequently designed to 
target this 150-cavity 22 . An open 150-loop has never been observed 
in X-ray crystallography experiments with group 2 NAs, especially in 
N2 containing a D147-H150 salt bridge, although simulation work 
suggested N2 may be able to adopt such a configuration in solution 
dynamics 21 . In this report, we demonstrate for the first time that the 
N2 150-loop can adopt a half-open conformation in its crystal 



structure. Furthermore we have carried out extensive molecular 
dynamics simulations in order to elucidate mechanisms of 150-loop 
motion. 

The free N2 crystal structure possesses a salt bridge between D147 
and HI 50, which is one of the key factors to determine the rigidity of 
its 150-loop. However, this rigid interaction can be broken under the 
condition of high concentrations oseltamivir carboxylate. With the 
exception of N5, which contains N147, all other types of NA possess 
G147. It can be inferred that the 150-loop should be more flexible 
with G147 in place of D147. Alternatively, the conformation of the 
150-loop in these NAs should be more easily influenced by extrinsic 
factors. It seems that both group 1 and group 2 NAs are able to adopt 
an open 150-cavity under certain conditions. From this perspective, 
the conformation of the 150-loop may not be a true group specific 
feature which is consistent with previous molecular dynamics simu- 
lation results 21 . This notion can also be supported by the 150-cavity 
deficiency in the atypical group 1 member, 2009 pandemic Nl 25 . 

It is worthy to note that an open 150-cavity in all types of NA has 
never been observed in the presence of zanamivir and it has been 
confirmed that zanamivir remains effective against oseltamivir-res- 
istant viruses. Based on analysis of the complex structures of N2 and 
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Table 3 | Data collection and refinement statistics 

Parameter FreeN2 N2-oseltamivir open form N2-oseltamivir closed form N2-G1 47-oseltamivir 
Data collecting 

Space group C222i C222j C222j C222i 

Unit cell dimensions (a, b,c) 114.83,39.47,140.05 114.91,139.39,140.24 115.43,139.11,140.00 114.33,139.31,140.05 

Unit cell dimensions (a, p, y) 90.00, 90.00, 90.00 90.00, 90.00, 90.00 90.00, 90.00, 90.00 90.00, 90.00, 90.00 

Resolution range (A)° 50.00-1.80(1.86-1.80) 50.00-1.60(1.65-1.60) 50.00-1.80(1.86-1.80) 50.00-2.20(2.28-2.20) 

R merge (%) b 10.2(54.7) 8.1(53.8) 9.7(49.2) 17.3(55.4) 

l/a 17.59(2.54) 19.1(2.29) 20.19(4.3) 11.17(3.53) 



Redundancy 
Completeness (%) 


6.2 (5.9) 
96.0 (98.1) 


4.9 (4.7) 
98.4 (93.9) 


6.6 (6.5) 
100 (100) 


7.2 (7.2) 
99.7(100) 


Refinement 


Resolution (A) 
Rwork (%)° 
Rfree (%) d 


37.45-1.80 
16.0 
17.8 


30.00-1.60 
13.7 
17.0 


38.75-1.80 
14.9 
16.9 


37.40-2.20 
19.9 
23.8 


RMSDs 


Bond lengths (A) 
Bond angles (") 


0.006 
1.199 


0.007 
1.231 


0.007 
1.219 


0.006 
1.245 


Ramachandran plot Quality 


Favored (%) 
Allowed (%) 
Outlier (%) 


96.5 
3.5 
0 


96.6 
3.4 
0 


96.5 
3.4 
0 


95.9 
4.1 
0 


"Value in parentheses refer to the highest 
b Rmerge = ihkl l-<l> | /Ehkll where 1 i 
°Rwork = S | Fobs - Fcalc | /SFobs where 
d Rfree is calculated as for Rcryst but usinc 


resolution shell of data, 
the intensity of unique relfection hkl and <l> is the average over symmetry-related 
Fobs and Fcalc are the observed and calculated structure factors, respectively. 
5% of reflections sequestered before refinement. 


observations of unique reflection hkl. 





inhibitors, we note that the 4-guanidino group of zanamivir provides 
stronger interactions with the NA active site to facilitate high affinity 
binding. This might compensate for lost interactions resulting from 
drug resistant mutations. A previous study illustrates that I223R (IR) 
and I223R/H275Y (IRHY) mutations of H1N1-2009 NA lead to 
reduce binding efficacy of oseltamivir carboxylate. The 150-loop of 
IRHY:osel shows a marked conformational change to the open form, 
resulting in a loss of the hydrogen bonds between the inhibitor and 
the 150-loop, contributing to resistance to oseltamivir 34 . Whearas the 
avidity of zanamivir is affected less by the IR and IRHY mutations as 
it is able to maintain hydrogen bonds to the 150-loop. Moreover, a 
R152K mutation in flu B results in oseltamivir and zanamivir resist- 
ance 30,31 demonstrating that the interactions between key 150-loop 
residues in the 150-loop and inhibitors play an important role in high 
avidity binding. Therefore, increasing the hydrogen bonds between 
NA inhibitors and residues in the 150-loop may help to overcome 
drug resistance. Additionally, it is difficult for the zanamivir guani- 
dino group to lose its charge due to its high pKa, and the positively 
charged guanidino group has a high affinity for the 150-loop. 

In conclusion, this study demonstrates the inherent flexibility of 
the 150-loop and the importance of inhibitor protonation state on 
influencing active site conformations. The open 150-cavity discov- 
ered in group 2 NA suggests that protein active site conformations 
can be greatly influenced by introduction of a ligand. Thus, both 
intrinsic protein dynamics and ligand induced protein dynamics 
should be considered in the future antiviral drug design. 

Methods 

Recombinant NA production. NA was prepared in baculovirus expression system as 
previous reported 1425 . cDNA encoding amino acid residues 83-469 of p57N2 (1957 
Asian pandemic H2N2 NA, Sequence GenBank accession No. M18648.1) were 
cloned into the baculovirus transfer vector pFastBacl (Invitrogen), with a GP67 signal 
peptide, a 6X histag, a tetrameric sequence and a thrombin cleavage site at the N 
terminus. A point mutation corresponding to G147 was introduced by site directed 
mutagenesis based on the p57N2 sequence. Recombinant baculovirus was prepared 
based on the manufacturer's protocol (Invitrogen). Sf9 cells were used to prepare 
recombinant baculovirus. 



Protein crystallization and inhibitor soaking assays. Protein concentrations were 
determined using a BCA protein assay according to the manufacturer's instructions 
(Pierce). Crystallization screens were initiated using commercial kits (Hampton 
Research). All crystallizations were performed using the hanging drop vapor diffusion 
technique. The N2 wild type and N2-G147 crystals were observed in PEGRx I 27 
(0.1 M BIS-TRIS propane pH9.0, 10% v/v Jeffamine ED-2001 pH7.0) and PEGRx 1 26 
(0.1 M Imidazole pH7.0, 20% v/v Jeffamine ED-2001 pH7.0), repectively. The N2 
crystals were incubated in mother liquor containing 20 raM and 40 uM oseltamivir 
carboxylate, and 20 mM Zanamivir at 291 K for 1 hour, 30 min and 1 hour, 
respectively. The N2-G147 crystal was incubated in mother liquor containing 20 mM 
oseltamivir carboxylate for 1 hour at 291 K. 

Data collection, structure determination, refinement and analysis. All crystals 
were soaked in cryoprotectant (reservoir solution supplemented with 20% glycerol), 
then flash-cooled and maintained at 100 K in a cryostream. Diffraction data for the 
free p57N2 native structure were collected at KEK beamline NE3A, while the data for 
N2 soaked with oseltamivr carboxylate and zanamivir were collected at SSRF 
beamline BL17U. Diffraction data were processed and scaled using HKL2000 35 . Data 
collection and processing statistics are summarized in Table 3. The structure of 
p57N2 was solved by molecular replacement using the A/TOKYO/3/1967 H2N2 N2 
molecule (PDB HVG) 36, as the search model with Phaser 37 in the CCP4 program 
suite 38 . The other datasets were determined by molecular replacement using the 
solved N2 structure. Extensive model building was performed by hand with COOT 39 , 
and restrained refinement was performed using REFMAC5 40 . Further rounds of 
refinement were performed using the PHENIX package 41 . The final structures 
displayed good stereochemistry (Table 3), as assessed by the program PROCHECK 42 . 
Structural figures were prepared with PyMOL (http://www.pymol.org/). 

Molecular dynamics simulations. Five simulation systems were built up and 
designated as OST-NH3-H150-Pos, OST-NH2-H150-Pos, OST-NH2-H150-Neu, 
ZMR-H150-Pos, and FREE-H150-Neu (Table 1). Atomic coordinates of systems 
OST-NH3-H150-Pos, OST-NH2-H150-Pos, and OST-NH2-H150-Neu were taken 
from the crystal structure of the N2- oseltamivir closed form complex; FREE-H150- 
Neu were taken from the crystal structure of the free N2 and ZMR-H150-Pos were 
taken from the crystal structure of the N2-zanamivir complex (PDB 3TIC) 27 . NA 
tetramers were used in each system. Protonation states of HI 50 and states of ligand 
were assigned as shown in Table 1. Water molecules and calcium ions in crystal 
structures were retained. 

RESP 43 charge of OST and ZMR were generated using antechamber 44 after 
calculation by Gaussian09 45 (HF/6-31 1G**). Force fields for OST and ZMR were 
generated based on general AMBER force field in amber package 43 . Then the system 
was setup using AmberTools. Amber ff99SB force field was applied to proteins 46 . The 
calcium was parameterized in the classical force field. TIP3P waters were added to 
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solvate each system with a boundary of 8 angstroms away from any protein or ligand 
atoms 47 . Counterions were added to neutralize each system. The system details can be 
found in the Supplementary information (Supplementary Table 1). Then the topo- 
logy of amber format was transformed to GROMACS readable format. 

Simulations were performed using the GROMACS 4.5.3 package 48 . To alleviate any 
steric clashes, the initial structures were minimized using the steepest descents 
method with a force tolerance of 10 kj mol -1 nm" 1 . Non-hydrogen protein atoms, 
main chain atoms and protein a-carbon atoms were progressively restrained with a 
force constant of 1000 kj mol -1 nm" 1 . After energy minimization, the systems were 
heated to 310 K with position restraints on non-water atoms. The temperature was 
controlled using the V-rescale method with a coupling time of t T = 0.1 ps and the 
pressure was maintained at 1 bar using Berendsen barostat with t p = 0.1 ps and a 
compressibility of 4.5 X 10~ 5 bar. Long range electrostatic interactions were calcu- 
lated using the particle-mesh Ewad (PME) algorithm 49,50 . The LINCS 51 algorithm was 
used to constrain hydrogen bond lengths. Further three 500 ps long equilibration 
runs were position restrained during the whole equilibration process. Finally, 200 ns 
MD production runs were performed for each of the five tetramer systems. The time 
step for integration was set to 2 fs. The coordinates of the systems were saved every 
10 ps. 

Population analysis based on 150-cavity. Conformations in each chain were 
extracted by an interval of 100 ps from the trajectories for 150-cavity volume 
calculation. The volume of 150-cavity was measured using the pocket volume 
measuring algorithm POVME 52 . The volume inclusion sphere that defines the 150 
cavity was a single sphere that encompassed the 150-cavity in the widely investigated 
group 1 neuraminidase (PDB 2HU0) 20 . Each conformation extracted was aligned to 
2HU0, and POVME algorithm neglected the volume occupied by atoms within the 
150-cavity sphere. 

Hydrogen bond statistics and RMSD calculation. VMD 53 was used to calculate the 
hydrogen bonds between two groups or residues. Hydrogen bond distance cutoff was 
set to 3.5 A and angle was set to 20 degrees. VMD was also used to calculate the RMSD 
of 150-loop in crystal structures. All a-carbon atoms were superimposed at first, and 
then the RMSD of a-carbon atoms in the 150-loop was calculated. Only Ca atoms 
were used to calculate RMSD values. 
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